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^ : ABSTRACT 
> 

We investigate the prospects for detection of lepton flavour violation in sparticle production and 
decays at a Linear Collider (LC), in models guided by neutrino oscillation data. We consider 
both slepton pair production and sleptons arising from the cascade decays of non-leptonic 
sparticles. We study the expected signals when lepton-flavour-violating (LFV) interactions are 
induced by renormalization effects in the Constrained Minimal Supersymmetric extension of the 
Standard Model (CMSSM), focusing on the subset of the supersymmetric parameter space that 
also leads to cosmologically interesting values of the relic neutralino LSP density Emphasis is 
given to the complementarity between the LC, which is sensitive to mixing in both the left and 
right slepton sectors, and the LHC, which is sensitive primarily to mixing in the right sector. 
We also emphasize the complementarity between searches for rare LFV processes at the LC 
and in low-energy experiments. 
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1 Introduction 



In recent years, a plethora of data from atmospheric pQ, solar [2] and long-baseline reactor [3] 
and accelerator [H [5] neutrino experiments have established the existence of neutrino masses 
and oscillations with near-maximal — v T and large v e — > mixing. A natural expectation 
in this context is that charged-lepton-flavour violation (LFV) should occur at some level. This 
may be enhanced sufficiently to become observable in a class of theories predicting new physics 
at the TeV scale accessible to colliders, particularly in supersymmetric theories. There are 
several sources of lepton flavour violation in such theories, which could have unacceptably large 
LFV if the soft supersymmetry-breaking masses of different sfermion flavours were not universal 
at some level. For this reason, it is often assumed that these masses are equal at the grand- 
unification scale, as in the Constrained Minimal Supersymmetric extension of the Standard 
Model (CMSSM). 

Even if the sfermion mass matrices are diagonal at the unification scale, as in the CMSSM, 
quantum corrections would modify this structure while running from the GUT scale to low 
energies. This effect is particularly interesting in see-saw models for neutrino masses, where 
the Dirac neutrino Yukawa couplings cannot be diagonalised simultaneously with the charged- 
lepton and slepton mass matrices [6j. Given the large mixing of the corresponding neutrino 
species, charged-lepton-flavour violation may occur at enhanced rates for sufficiently small 
soft supersymmetry-breaking masses, giving rise to observable signals such as fi — > ey, \i — e 
conversion, r — )■ fij and r — > EJ E] . 

Other charged-lepton-flavour violating possibilities that have been considered include slepton 
pair production at a Linear Collider (LC) PHI ELTJ 021 E2 EH], and also signals at the LHC 
[To"! ITo^ [T7] . particularly in % 2 — > X + e ± /i =F X2 — > X + Ai ± r T decays (here x is the lightest neu- 
tralino, assumed to be the lightest supersymmetric particle (LSP), and \2 is the second-lightest 
neutralino). These decays potentially provide search prospects that are complementary to di- 
rect searches for the flavour-violating decays of charged leptons. However, most of the previous 
studies of LFV at a LC use a low- mass spectrum (focusing on P cm < 400 GeV or y/s < 800 GeV) 
and the points chosen for study do not always satisfy all the relevant phenomenological and 
experimental bounds (e.g., [IE]), which are now being improved by the LHC [191 ED]- Moreover, 
previous analyses have not always taken into account the constraints on the cosmological relic 
LSP density imposed by WMAP and other experiments [21 J . 

In the current paper, we revisit the various detection channels for lepton flavour violation at 
a Linear Collider, studying the complementarity between different processes. In all cases, we 
focus on regions that satisfy not only all phenomenological constraints, but also the cosmological 
relic density considerations. We extend previous results to somewhat heavier sparticle spectra, 
and make comparisons with the expectations for experimental sensitivity at the LHC. 

Our paper is structured as follows. In Section 2, we briefly review the theoretical framework 
and sources of slepton mixing and summarize the possibilities for observable signatures in 
slepton production at a LC. In Section 3, we discuss the relevant supersymmetric parameter 
space, and pick representative points motivated by both phenomenological and cosmological 
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considerations. In Section 4, we look at the expected slepton mixing parameters, with emphasis 
on the region where LC experiments could give information additional to that obtainable from 
low-energy LFV decays. In Section 5, we discuss in more detail the cross sections, comparing 
with the expected signals at the LHC, where appropriate. Finally, in Section 6, we summarize 
our conclusions and outlook. 



2 LFV in Slepton Pair Production 

2.1 CMSSM with LFV 

In the unrotated charged-lepton flavour basis li = (e^, £ll, ?l, jl R , f R ), the charged slepton 
mass matrix is: 

.,2 / M 2 M 2 \ , 1N 



where 



* 1^ M RR 



M 2 LL = m\m e + Ml - ~(2m 2 w -m 2 z ) cos 2(3 I, 



M 2 RR = m\m t + M R - (m| - m 2 w ) cos 2/3 I, 
m lr = (A e - fx tan (3) me, 

M 2 RL = (M 2 R y. (2) 

Here we parametrize trilinear soft supersymmetry-breaking terms as Afj ■ (\ e )ij, where the A e 
are the respective Yukawa couplings. For universal soft terms at some high input scale, one has 

M\ = M 2 R = m 2 J, A% = A ■ %, (3) 

whereas flavour-mixing entries may be parametrized by: 

4x = {M 2 xx T/{M 2 xx f (X = L,R). (4) 

The correspondence between the mixing parameters 5 XX and the flavour mixing parameters 
used in the phenomenological study of slepton production cross sections is easily derived. For 
example, defining the splitting of the third-generation soft supersymmetry-breaking mass soft 
terms as might be generated by the renormalization-group equations (RGEs) by 

K xx = l- (M 2 XX ) 33 /(M 2 XX ) 22 , (5) 

the slepton mixing angle 823 and the splitting of the second- and third-generation mass eigen- 
states Arfi23 can be written as: 

tan{26 23 ) = 25 XX /K XX , (6) 
m sin(29 23 ) 

The evaluation of the LFV observables is done by performing the diagonalization of the slepton 
mass matrices (see [8], for instance), inserting the full rotation matrices in the lepton-slepton- 
gaugino vertices and summing over all the mass eigenstates of the exchanged particles. 



3 





Figure 1: Feynman diagrams for e + e — > — > ^t^ a x\x\- The arrows on scalar lines indicate the 
flow of lepton number. Similar diagrams - appropriately modified - exist for charginos. 



2.2 LFV Cross Sections 

As already discussed in the introduction, charged-lepton flavour violation at a LC may oc- 
cur either directly in slepton pair production or indirectly via slepton production in cascade 
decays [TO] . Processes leading to lepton production in the decays of a pair of sleptons include 



\fq _> rVx?Xi, 



e+e ViV c - r^fi^XiXi , 



(8) 



for which representative Feynman diagrams are shown in Fig. 1. Slepton production may also 
result from the the cascade decays of the heavier gauginos, e.g., via the processes 



e + e" ->• xtXi -> ^^XiXi- 



(9) 



These amplitudes differ in the interference terms, and read as follows: 



Mlf = Mf r 



q 2 — mf + irhi 



-T in M+- 



»#casc \ * n/fcascrp 

1V1 aB — 1VI P J " 



q 2 — m 2 + irriiTi 



Ti ia D p 2 -m 2 i +im l T i ! > 

-T*nM 



(10) 

(11) 



for pair production and cascade slepton production, respectively, where Mp and Md are the 
production and decay amplitudes for sleptons in the absence of LFV, and T ia parametrizes the 
lepton- flavour- mixing matrix element. 

Complete expressions for the respective cross sections are given in [TO], and used in our work. 
However, intuition may be developed by first looking at the simplified case including only 2-3 
slepton mixing. In this case, the only diagrams contributing to the cross sections are those 
mediated by 7 and Z° exchange. Even simpler formulae can be obtained in certain limits, 
such as the cases of narrow widths and small mass differences between the sleptons of different 
generations. For example, if Am^ <Cm = \{ m 2 + anc l m ^ij — (m^i + rhjTj)/2 <C m 2 , in 
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the presence of 2-3 mixing the cross sections for pair and cascade slepton production are given 
by [22]: 

^ pair = X2 3 (3 - 4 X23 ) sin 2 2 6 23 a(f f t t)Br{£+ -> t X)Br(t -^TY), (12) 
<f = X23 sin 2 2^ 2 3 a(f f^tX t)Br(t -> f Y), (13) 

where a(f f — » £ + X£~), a(f f — > £ + fr) and Brii.^ — » ^X) are the cross sections and 
branching ratios in the absence of flavour violation, and LFV enters via X23 = #23/2(1 + x| 3 ) 
where 223 = Arri23/r23, and sin 2 2^23 parametrizes the slepton mixing angle. In the limit 
#23 3> 1, interference can be neglected and the results are even further simplified. However, in 
the case that the interference term dominates, it suppresses LFV significantly. 

The exact experimental signatures of the above effects depend on the decay chains, as discussed 
in previous analyses [Uj. Channels involving the lightest charginos are potentially interesting, 
particularly those with hadronic chargino decays, since the lightest neutralino is 'seen' as missing 
energy. In general, the dominant chargino decay is xt ~~ whereas channels such as xt ~^ 
q u qdX° are relevant when m~+ < or, in other words, when the xf is the NLSP. Channels 
leading to hadronic xt decays [11] are mediated only by left sleptons, whilst <r(e + e~ — > r ± y u T + 
2x°) may be mediated by either right or left sleptons. However, points with m~+ < are 
disfavoured in the framework of the CMSSM. For instance, the previously proposed point with 
tan/3 = 3, m = 100 GeV, m 1/2 = 200 GeV [12] was excluded by the LEP m h bound [IB] 
and is now incompatible also with LHC limits [19]. This could perhaps be achieved in the 
focus-point region, where rayi < 1 TeV and m > several TeV, but slepton masses in this 
region are very heavy and hence LFV processes are more suppressed. 

The channel e + e _ — > t + 2%° was studied in [13], where it was found that a signal of 1 
fb could be distinguished from the background. Only slepton mixing in the left sector arising 
from see-saw neutrino mass models was included in [TJ] , and the limit from r — > ^7 was at the 
time less restrictive. As shown in [T7], reasonable rates of LFV processes at the LHC in the 
slepton sector are possible while preserving the bounds on r — > /ry if one assumes mixing also 
among the right sleptons, as may occur in non-minimal GUT theories. 

3 Supersymmetric Parameter Space 

From the above discussion, it is clear that the dependence on the supersymmetric parameter 
space is particularly crucial, since the chains of cascade decays that will dominate, as well as 
the relative strengths of the signals and backgrounds, depend on the mass hierarchies of the 
superpartners. 

We start with some preliminary considerations of the relevant channels within the CMSSM pa- 
rameter space. In previous work [TTJ [23] a region that is promising from the cosmological point 
of view was discussed, in which sizable f — x coannihilations or direct-channel H/A resonances 
lead to values of £l x h 2 consistent with WMAP. As a representative of the coannihilation case, 
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we chose tan/3 = 35, a value that leads to a relatively heavy sparticle spectrum compatible 
with WMAP. At tan /3 = 45 we found points in the resonance funnel, which have larger scalar 
masses since m can take higher values than in the coannihilation area. We did not consider 
the focus-point region of the CMSSM because the superheavy slepton masses there lead to 
reduced LFV effects. We note that these regions are compatible with the LHC constraints on 
supersymmetry, which indicate that m\ii > 450 GeV and favour tan/3 > 15 |20j. 
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Figure 2: Contours of the cross sections a(e + e~ — > = 1 fb in the (mi/2,mo) planes for tan f3 = 

35 (left) and 45 (right) at three different energies, namely yfs = 500 GeV (red thin lines), y/s = 
1000 GeV (black thick lines) and yfs = 2000 GeV (blue thin lines). The upper and lower near- 
horizontal lines correspond to the following slepton pair production processes, respectively: f^f^ (solid) 
and v T u T (dash). From left to right, the near-vertical lines correspond to a{e + e- -> X \X°2) = 1 fb (dot- 
dash) and o-(e + e~ — > X1X2) = 1 ft* (double-dot-dash). In regions above the line where = m„ T 
(purple dots), chargino-mediated processes are allowed. Representative contours of BR(b — > sj) are 
indicated, and the regions with LSP relic density in the range preferred by WMAP are shaded green. 
The benchmark points selected for further study are indicated by crosses. 

In Fig. |2] we present contours of the production cross sections for sparticles that are relevant 
for our study of LFV effects within the CMSSM, using 1 fb as a reference cross-section value. 
We indicate in green the constraint on the parameter space derived from the WMAP-favoured 
CMSSM parameter space (0.094 < Vth 2 < 0.128 at 3 - a) using MicrOMEGAS [23], and in grey 
the regions excluded by the constraints that BR(6 — > 57) > 2.5 ■ 10~ 4 and > 111 GeV 
(we include an 3 GeV allowance for the uncertainty in the theoretical calculation of m^). We 
present contours for the production of fj rf calculated using CalcHEP [25J: the cross sections 
for producing other right sleptons would have similar contours, though at somewhat smaller 
values of (mi/ 2 ,m ) because they are heavier than the f\. The left slepton production cross 
sections are smaller, in general, and here we draw the a(p T v T ) = 1 fb contour (the cross sections 
for all slepton pair-production processes are compared in Fig. |3]). The relevant channels for 
LFV in cascade decays following pair production are also displayed in Fig. |2j along with the 
line m~+ = m^, above which chargino-mediated LFV processes are allowed. 
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Focusing initially on the lines corresponding to \fs = 500 GeV, we see that on-shell production 
of left sleptons would occur mostly in the region excluded by b — > 57. At this energy, therefore, 
we expect right slepton production to be more important. However, production of charged 
slepton pairs in the cascade decays of neutralinos is also possible, whilst chargino decays are 
disfavored. At \fs = 1000 GeV, we see a significant part of the parameter space with left slepton 
pair production (below the dashed line) and with both left and right slepton production (below 
the solid line). The production of pairs of charged sleptons in chargino cascade decay is not 
possible in most of the WMAP favored area, since m-+ < m^. Finally, at \fs = 2000 GeV all 
the decay channels are present in the region of parameter space favoured by WMAP. 

Based on this discussion, two benchmark points in the allowed regions for each of the choices 
tan/3 = 35,45 have been selected for subsequent study, namely: 



(a35) 


tan(3 


= 35 


m 


= 255 GeV, 


M 1/2 


= 610 GeV, 


(635) 


tanj3 


= 35 


m 


= 345 GeV, 


Mi/ 2 


= 900 GeV, 


(a45) 


tanfi 


= 45 


m 


= 440 GeV, 


Mi/2 


= 550 GeV, 


(645) 


tan(3 


= 45 


m 


= 730 GeV, 


Mi/ 2 


= 900 GeV 



with A = in each case, resulting in the particle spectra shown in Table 1. We note that 
all these points lie beyond the regions of CMSSM parameter space excluded by 2010 LHC 
data [HI [20]. In general, these points lead to more promising detection prospects at a LC at 
high energies: at \/s = 500 GeV, only production through right sleptons is allowed. 



Point 


m 


Mi /2 


tan /3 


M 4 




M x1 




M fl 


M f2 




m !l 


M, T 


M 5e 


a35 


255 


610 


35 


480 


741 


255 


480 


261 


477 


342 


478 


454 


472 


b35 


345 


900 


35 


720 


1034 


383 


720 


385 


670 


478 


683 


655 


679 


a45 


440 


550 


45 


432 


671 


229 


433 


338 


543 


485 


571 


519 


566 


b45 


730 


900 


45 


725 


1015 


385 


725 


590 


872 


801 


936 


858 


933 



Table 1: Relevant sparticle masses (in GeV) for the selected benchmark points. The selected points 
all have Aq = 0. 



We show in Fig. [3] the dependence of the flavour- conserving (FC) cross sections on the y/s of the 
collision for the benchmark points ( fl5l) with tan/3 = 35: the cross sections for the benchmarks 
with tan /3 = 45 are very similar. We can already get an idea of the LFV cross sections from the 
magnitudes of these FC cross sections, by looking at the approximate formulas ( TL21) and ( ITBl . 
which assume that only 2-3 mixing is present. These provide an initial picture of the order 
of magnitude of LFV (up to 40% error), simply by convoluting the slepton flavour- violating 
decays with the FC results. The latter can be extracted automatically, e.g., from CalcHEP |25j, 
allowing one to explore quickly the relevant parameter space. If the mixing was indeed limited 
to the 2-3 sector, the accuracy would be larger, within 10% of the exact calculation discussed at 
a later stage. Taking into account the full mixing (as shown below) leads to additional channels 
and an enhancement in the cross sections cx(e + e~ — > £i£j), also for the case of mixing in the LL 
sector, despite the fact that mixing with the first generation is very constrained by the bound 
BR(/x ->■ e7)< 1.1 • 10" 11 |2E]. 
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Figure 3: Energy dependence of the cross sections for sparticle production in the channels and 
(TJ|). The black thick (blue thin) lines correspond to point a35 (b35) in M5\) . The upper left (right) 
panels display left (right) charged slepton production, and the lower left (right) panels display sneutrino 
(gaugino) production. 

The cross sections e + e~ — > £^£f, % — L, 1, R, 2, which are relevant for calculating e + e~ — > 
Ijtj — > T ± /i =F XiXi) are displayed on the upper panels of Fig. [3j We see that the production 
rates for left and right sleptons are comparable. 

In the case of e + e~ — > z>jZ>| — > T ± /i T x^x ] ~, Fig. |3] indicates that the production rate for sneu- 
trinos can be larger than that of charged sleptons. However, i) even if e + e~ — > u^Uj is large, 
the decay Vi — > l\i is allowed only in a restricted area of the parameter space, as we can see 
in Fig. [21 and ii) the advantage of producing the Xi is n °t so important in this case, because 
the decay xt — > x + QuQd is not allowed in the range of tan/3 we favour (we recall that in the 
CMSSM the NLSP is usually a stau). 
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Finally, we comment on the cascade slepton production processes of These cases are not 
so promising within the CMSSM, but the analysis is simple, using the flavour-conserving cross 
sections for charginos and neutralinos and multiplying by the LFV branching ratio. The LFV 
signal arises from the decay of the X2 or xt > which are produced at lower rates than the sleptons. 
The LFV decay of the X2 has been already studied in the context of the LHC in [17]. Here, in 
agreement with [TTJ, we find that mixing in the RR sector may lead to larger signals. However, 
the resulting cross sections are smaller than those mediated by sleptons. Also, the channel 
with xt is n °t so relevant, not only because the production of charginos is lower than that of 
neutralinos (as can be seen in Figs. |3]), but also because it is only sensitive to LL mixing in 
points where m-+ < m„ T . 



4 Opportunities for LFV observation at the LC 

LFV effects are severely constrained by decays of the form BR(£j — > £jj), which limit signifi- 
cantly the allowed parameter space, with respect to both the mixing terms 5ll,rr and the values 
of the soft supersymmetry-breaking masses. These limit the range of parameters in which we 
can get observable cross sections without violating the LFV bounds. In general, we may expect 
that LFV decays are progressively suppressed for a heavier sparticle spectrum at fixed values of 
8ll,rr, and hence that one may get observable effects only for larger values of 5ll,rr- However, 
since the LC gives direct access to energies that are relatively high, the suppression induced by 
a heavier spectrum may be smaller than those in rare LFV decays and /1 — e conversion. It is 
natural, therefore, to ask the following questions, (i) For a heavy sparticle spectrum, how small 
are the 5ll,rr that are accessible before the cross sections become too small to be observable? 
(ii) For a light sparticle spectrum, what are the minimum values of Sll^r that can be probed? 

These issues are addressed in Fig. HI where we display values of the flavour- violating parameters 
5rr and 5ll, showing lines of constant a(e + e~ — > r ± /i T + 2%°) for y/s = 2000 GeV0. Using 
the bounds on r — )■ /ry and r — > /ry, we display the allowed ranges in ((6^)23, (5ll)i3) planes 
(upper panels), and the corresponding mixing parameters in the RR sector (lower panels). In 
this exploration, the values of rr)i2 were chosen to be proportional to those of (5ll, k_r)i3, 
so as to reduce the number of parameters. We have checked that the dependence of our results 
on these parameters are small, provided they are taken in a range such that the bound on 
BR(jU — > ej) is preserved. In addition, we note that the ratio (Sll,rr)i3 ((<*> ll,rr) 12) has a very 
small influence in r — > ^ (t —> ej). 

As we see in Fig. HJ the constraints from the LFV branching ratios result in horizontal and 
vertical lines in the 5ll and 5rr planes, corresponding to a box of allowed parameters that 
scales with the masses in the spectrum. On the other hand, the lines of constant LFV cross 
sections depend strongly on the kinematical properties of the chosen point as well as on the 
size of the mixing. In particular, we note the following: 

1 We choose y/s = 2000 GcV in Fig. 2] because, had we chosen y/s = 1000 GcV, only the RR channels would 
lead to on-shell sleptons. Even for = 2000 GeV, we can see in Fig. [5] that point &45 predicts very low cross 
sections for the LL case, due to the fact that the left sleptons are only marginally accessible, kinematically. 
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Figure 4: Constraints on the magnitudes of the mixing parameters and possible LFV effects for points 
with tan (3 = 35 (left panels) and tan/3 = 45 (right panels). The shaded areas are those allowed by 
current limits on BR(t — > cy) (dot-dash line) and BR(t — > /j,j) (dash line) using a35, a45 as reference 
points (thick lines bounding the solid shaded areas) and , 635, 645 (thin blue lines bounding the ruled 
shaded areas). The solid lines are contours of o~(e + e~ — > t (jF + 2^°) in fb for y/s = 2000 GeV. 

• The shapes of the lines are quite different for each benchmark point, as was to be be expected 
from the variations of the cross sections between the different benchmark points in the WM AP 
favoured area. For instance, for LL mixing, the differences between the predictions for the two 
points at tan (3 = 45 (upper right panel in Fig. HJ) can be understood from the behaviours of 
the cross sections shown in Fig. |5j 

• Except for LL mixing at tan /3 = 45, we obtain LFV cross sections in the fb range at points 
of the parameter space where LFV decay branching ratios are below the present experimental 
bounds. This reflects the observation that passing to a heavier spectrum increases significantly 
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Figure 5: Values of o~(e + e~ — > t^tj — > r ± // =F + 2x°) vs \fs for the benchmark points with tan/3 = 35 
(left) and tan/3 = 45 (right): solid (dashed) lines correspond to a35 and a45 ( 635 and 645 ). Along 
the thick lines we assume the RR mixing of U6\) . and along the thin lines the LL mixing of H5\). 

the allowed size of the mixing parameters, due to the rapid decrease of the LFV branching 
ratios for heavy sparticle spectra. 

• It is interesting to remark that at the LC mixing in the LL sector is observable without 
exceeding the bounds on £ { ij'j |26j, due to the presence of additional channels as compared 
to the LHC In the case of the LHC [17], mixing in the RR sector (which is less constrained 
but would require a departure from the CMSSM) had turned out to be the most promising 
avenue for discovering flavour violation. This is a clear distinction between the two colliders, 
and implies an advantage for the LC, in that its searches can be more directly connected with 
neutrino mass and mixing parameters. 

The above results motivate a particular choice of parameters for a more detailed discussion of 
LFV cross sections. We assume the same values of the mixing parameters for both choices of 
tan/3, namely: 

= 0.02, (5 LL ) 23 = 0.02, (15) 

OMia = 0.04, (<W)23 = 0.15. (16) 

For completeness, we also introduce a small mixing between the first and second generation: 
{8 ll/rr) 12 = 0.2 • (Sll/rr)i3- The LFV decays into r — fx pairs are not heavily dependent on 
this parameter while, with this choice, the bound on BR(/i — > ey) does not over-constrain the 
parameter space. 

LFV in the fi — e channel is generally suppressed in our scheme, unless 5lr mixing is introduced. 
In our analysis, however, we do not include such a mixing, in order to be able to discuss 
separately the LFV arising from each chiral sector, thereby avoiding the additional model 
dependence associated with choosing 5lr ^ 0. In the case of large RR mixing, non-trivial 
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interference of the RR sector with the LL sector would arise through LR mixing. The // — > e^y 
bound, which could almost be ignored in selecting RR and LL mixing, would now play a very 
important role, forcing us to either a very small part of the parameter space or to very small 
values of (Slr)- To reach the bound on radiative decays with LR contributions alone we would 
need: (5 LR ) 12 ~ 10- 3 (<^)i3,i3 ~ 10~ 2 [27]. 

We also followed the conservative approach of keeping A = 0, avoiding the introduction of other 
new LFV parameters. The main consequence for our analysis of allowing A ^ would, since 
the first and second generations have small Yukawa couplings, be to allow the chiral composition 
of the staus to vary. However, additional LFV could be introduced by assuming a non-trivial 
mixing matrix for the A terms. In this case, chiral mixing would occur for all generations, 
restricting among other effects the size of the RR mixing. The decay \i — > cy is the most 
sensitive to LFV A terms, and establishing a correlation between /i and r decays without having 
an specific model would be non-trivial. Nevertheless, since the bound on BR(/i — > ej) imposes 
a very severe restriction on on (M LR )i 2 , even assumptions such as (M LR ) 23 ~ (M LR ) ri ~ 
100(M LR )i 2 would not alter our main conclusions. 

5 Study of LFV at the LC 

The next step is to understand the dependence of the results on the available y/s. Clearly, at 
the LC, higher values of yfs provide accessible cross sections for heavier sparticle spectra, a 
potential advantage of the LC over the LHC, where it is relatively difficult to reach large values 
of the scalar masses. The variation of <r(e + e~ — > r ± /i =F + 2%°) with yfs is displayed in Fig. [5j 
where it is shown how, as yfs = 2P cm increases, it becomes possible to produce on-shell sleptons 
with larger masses. We should underline once more that the effects due to the kinematics as 
relevant as those from the LFV mixing parameters. The shapes of the lines can be understood 
as follows. 

• The cross sections for slepton pair production show interesting dependence on yfs, since addi- 
tional channels arise from slepton mixing. We see in Fig. [5] steps in the left slepton production 
cross sections, which appear as the sleptons of the first two generations become kinematically 
accessible. The steps corresponding to lighter generations are less pronounced in the case of 
RR mixing than in the LL case. 

• The overall cross sections for slepton production in in cascade decays are independent of the 
LFV parameters, and LFV arises in the decays of the gauginos. Figs. [6] and [7J display the cross 
sections for the pair production of x 2 Xi or X2X1 multiplied by the LFV branching ratios for 
the decays of x 2 or X2 1 computed as in [T7j. This explains the simple shapes in this case. 

• Distinguishing between pair and cascade slepton production is quite crucial, particularly since, 
as already mentioned, it is a fundamental difference between the LC and the LHC. Comparing 
the pair and cascade slepton production cross sections in Figs. and M (where we chose similar 
scales to facilitate comparisons) we see the following. In slepton pair production, both LL 
and RR mixing effects are of the same order of magnitude. On the contrary, in the case of 
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Figure 6: Same as in Fig. but for cascade slepton production in the channel a(e + e — > X1X2 
t ± /^ + 2 X °). 




Figure 7: Left panel: Same as in Fig. [5|, but for cascade slepton production in the channel a{e + e~ — > 
X1X2 ~> t^/x^ + 2%°). Right panel: Same as in Fig. O but for pair production with sneutrino 
mediation a(e + e~ — > UiUj —> r =t ^ =F + XiXi)- F° r tan/3 = 35 these channels are not allowed because 
mo < . 

cascade slepton production, the channel with LL mixing has a cross section at least one order 
of magnitude less than that with RR mixing, offering a distinction between the two mixing 
scenarios. 

• Channels with charginos in the final states could be promising when m~+ < m^,, i.e., when 
the xt is the NLSP. However, as already remarked, this is not the case of the CMSSM. In 
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general, the dominant chargino decay is xt —> lv, thus we do not obtain additional information 
beyond that from simple slepton production and decay. This is true also for channels with the 
heavier neutralinos. 

• Cascade slepton production mediated by charginos is possible only with mixing in the LL 
sector. Fig. [TJ (left panel) indicates that this channel could be relevant at y/s larger than in 
the neutralino mediated channels, and the cross sections would be of comparable magnitude 
in the LL mixing case. On-shell production requires irty > m-+, which is not the case for the 
sample points with tan (5 = 35, as they lie in the coannihilation region of the WMAP allowed 
area. However, for tan (3 = 45 the WMAP region includes points with resonant annihilation of 
neutralinos where the condition > m~+ is satisfied. 

• The expectations for sneutrino pair production and subsequent chargino production are shown 
in the right panel of Fig. [7J Relaxing the assumptions of the CMSSM can change the kinematics 
of cr(e + e _ — » X1X2), as well as the mass ratio of xt versus v that determines the decay 

Cross sections for e + e~ — » r ± /i =F + 2%° above 1 fb can be reached in the areas of parameters 
shown in Fig. [8] for energies below those indicated. The large areas between the thick solid 
red, black and blue lines, on the one hand, and the shaded regions that are excluded by the 
indicated present experimental constraints, on the other hand, demonstrate that there are 
ample opportunities for LFV discovery and measurement at the LC. These opportunities are 
exemplified by the benchmark points marked by crosses in Fig. |SJ 

In Fig. [9] we present the maximum values of the cross sections in the areas allowed by all the 
current constraints. The shaded areas show the possible ranges of <r(e + e~ — » r ± /i T + 2%°), 
and we note that points along the WMAP strips (solid lines) generally have high values of the 
cross sections. The dashed lines show the possible values of a(e + e~ — > r ± e T + 2%°) along the 
WMAP strips, and we see that these cross sections may be of the same order of magnitude as 
those for /i — r pairs H 

In Fig. [10] we display the expectations for the production of LFV r — [A pairs in the cascade 
decays (J5J), which are suppressed compared to the pair production process (jHJ). The production 
of LFV pairs from X2 decays is analogous to the corresponding process at the LHC [T7J We 
observe that phenomenologically interesting values on the LL sector induce very low cross 
sections, whereas mixing in the RR sector may be more observable, as it enhances the expected 
rates. 



6 Conclusions 

Motivated by the challenging possibility of observable charged lepton flavour violation, which 
would provide a new window on flavour physics that illuminates novel aspects of supersymmetry, 

2 We do not display values for a(e + e~ — > /i ± e =F + 2%°), because this would require a more precise analysis of 
mixings in the LR sector. 
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Figure 8: The solid red (black) (blue) lines are contours where a(e + e~ — > i /x T + 2x°)= 1 fb at 
yfs = 500, 1000, 2000 GeV, assuming the left mixing parameters [15\) in the upper panels and the 
right mixing U6\) in the lower panels. The left panels are for for tan /3 = 35 and the right panels 
for tan/3 = 45, assuming Aq = in both cases. Each panel also shows the areas excluded by current 
bounds on BR(/i — > ej) (thin solid line), BR(t — > e^f) (thin dash line), BR(t — > ^) (thin dot-dash 
line). The areas excluded by BR(b — > sj) and the LEP Biggs search are also displayed, and the green 
area denotes the WMAP favored region. We see that there are ample opportunities for LEV discovery 
and measurement at the LC, and the benchmark points chosen for further studies are indicated by 
crosses. 



we have explored possible signatures at the Linear Collider (LC), combining phenomenological 
constraints with the cosmological considerations implied by WMAP. We have found that both 
direct slepton pair production and sleptons produced in cascade decays may provide interesting 
signals in the cosmologically-favoured region of the supersymmetric parametric space. More- 
over, the LC could provide additional insights beyond those obtainable from the LHC, by virtue 
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Figure 9: The shaded areas show the possible ranges of a(e + e~ — > r ^ + 2%°) for tuq < 1000 GeV, 
tan/3 = 35 (left panels) and tan/3 = 45 (right panels), assuming Aq = in both cases, with ^/s 
fixed to 500 GeV (grey), 1000 GeV (green) and 2000 GeV (orange). In the upper panels, we assume 
the left mixing parameters ( fi5j) . whereas in the lower panels we assume the right mixing (16]) . The 
solid lines present the possible values for models along the center of the WMAP strips, with the lines 
corresponding to ^/s=500, 1000, 2000 GeV being progressively thicker. The corresponding predictions 
for a(e + e~ — > r ± e T + 2%°) are shown by the dashed lines. 

of its greater kinematic range for slepton production and its sensitivity to LL mixing as well 
as RR mixing. 

Within this framework, we found the following. 

• The LC enhances significantly the prospects of detecting LFV for heavy sparticle spectra, 
where favour-violating rare decays and conversions are significantly suppressed. This allows 
probing an entirely different range of the flavour-violating parameters. 
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Figure 10: Predictions from cascade decays (TJ|) along the WMAP strips for yfs= 1000 GeV (thick 
black lines) and 2000 GeV (thin blue lines). These channels are not accessible at ^/s=500 GeV. The 
solid (dashed) lines correspond to mixing in the RR (LL) sectors. 

• Unlike the LHC, at the LC we can have significant lepton flavour violation within the CMSSM 
via LL mixing. 

• The decay a(e + e~ — > r ± fi T +2x°) provides an optimal search channel in this respect. Channels 
where charginos play a significant role may require a departure towards theories where gaugino 
or Higgs unification is broken, in order to respect the conditions imposed by cosmology. 

• Comparing direct slepton pair production with indirect slepton production through the cas- 
cade decays of heavier sparticles, we find that in the latter case, the cross sections induced by 
LL mixing are approximately one order of magnitude lower than in the RR mixing case. This 
is consistent with the results previously obtained for the LHC, and is due to the presence of 
additional channels in cascade pair production. These lead to an interesting scaling of the cross 
sections with the available energy, which is sensitive to the 3x3 slepton mixing parameters. 

It would be very interesting to further investigate the following issues. 

• The cross sections expected in non-minimal extensions of the theory might not only enhance 
channels that in the current scheme are more suppressed, but could also enable a comparison 
of the allowed range of mixing parameters in different models. 

• The expectations for models of massive neutrinos, in which quantum corrections provide a 
significant source of LFV in the LL channel, provide a potential link between LC observables 
and neutrino mass and mixing parameters. 

• Detailed simulations of signals and backgrounds are desirable. 

Overall, it seems that the LC provides an optimal environment for the study of LFV, whereas 
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the LHC is limited to specific channels that have significant backgrounds. The fact that the LC 
opens up additional possibilities may prove significant for making the link between observable 
cross sections and flavour model building. 
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^ : ABSTRACT 
> 

We investigate the prospects for detection of lepton flavour violation in sparticle production and 
decays at a Linear Collider (LC), in models guided by neutrino oscillation data. We consider 
both slepton pair production and sleptons arising from the cascade decays of non-leptonic 
sparticles. We study the expected signals when lepton-flavour-violating (LFV) interactions are 
induced by renormalization effects in the Constrained Minimal Supersymmetric extension of the 
Standard Model (CMSSM), focusing on the subset of the supersymmetric parameter space that 
also leads to cosmologically interesting values of the relic neutralino LSP density Emphasis is 
given to the complementarity between the LC, which is sensitive to mixing in both the left and 
right slepton sectors, and the LHC, which is sensitive primarily to mixing in the right sector. 
We also emphasize the complementarity between searches for rare LFV processes at the LC 
and in low-energy experiments. 
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1 Introduction 



In recent years, a plethora of data from atmospheric pQ, solar [2] and long-baseline reactor [3] 
and accelerator [H [5] neutrino experiments have established the existence of neutrino masses 
and oscillations with near-maximal — v T and large v e — > mixing. A natural expectation 
in this context is that charged-lepton-flavour violation (LFV) should occur at some level. This 
may be enhanced sufficiently to become observable in a class of theories predicting new physics 
at the TeV scale accessible to colliders, particularly in supersymmetric theories. There are 
several sources of lepton flavour violation in such theories, which could have unacceptably large 
LFV if the soft supersymmetry-breaking masses of different sfermion flavours were not universal 
at some level. For this reason, it is often assumed that these masses are equal at the grand- 
unification scale, as in the Constrained Minimal Supersymmetric extension of the Standard 
Model (CMSSM). 

Even if the sfermion mass matrices are diagonal at the unification scale, as in the CMSSM, 
quantum corrections would modify this structure while running from the GUT scale to low 
energies. This effect is particularly interesting in see-saw models for neutrino masses, where 
the Dirac neutrino Yukawa couplings cannot be diagonalised simultaneously with the charged- 
lepton and slepton mass matrices [6j. Given the large mixing of the corresponding neutrino 
species, charged-lepton-flavour violation may occur at enhanced rates for sufficiently small 
soft supersymmetry-breaking masses, giving rise to observable signals such as fi — > ey, \i — e 
conversion, r — )■ fij and r — > EJ E] . 

Other charged-lepton-flavour violating possibilities that have been considered include slepton 
pair production at a Linear Collider (LC) PHI ELTJ 021 E2 EH], and also signals at the LHC 
[To"! ITo^ [T7] . particularly in x° — >* X + e ± /i =F ~* X + /W ± r T decays (here x is the lightest neu- 
tralino, assumed to be the lightest supersymmetric particle (LSP), and \2 * s the second-lightest 
neutralino). These decays potentially provide search prospects that are complementary to di- 
rect searches for the flavour-violating decays of charged leptons. However, most of the previous 
studies of LFV at a LC use a low- mass spectrum (focusing on P cm < 400 GeV or y/s < 800 GeV) 
and the points chosen for study do not always satisfy all the relevant phenomenological and 
experimental bounds (e.g., [IE]), which are now being improved by the LHC [191 ED]- Moreover, 
previous analyses have not always taken into account the constraints on the cosmological relic 
LSP density imposed by WMAP and other experiments [21 J . 

In the current paper, we revisit the various detection channels for lepton flavour violation at 
a Linear Collider, studying the complementarity between different processes. In all cases, we 
focus on regions that satisfy not only all phenomenological constraints, but also the cosmological 
relic density considerations. We extend previous results to somewhat heavier sparticle spectra, 
and make comparisons with the expectations for experimental sensitivity at the LHC. 

Our paper is structured as follows. In Section 2, we briefly review the theoretical framework 
and sources of slepton mixing and summarize the possibilities for observable signatures in 
slepton production at a LC. In Section 3, we discuss the relevant supersymmetric parameter 
space, and pick representative points motivated by both phenomenological and cosmological 



2 



considerations. In Section 4, we look at the expected slepton mixing parameters, with emphasis 
on the region where LC experiments could give information additional to that obtainable from 
low-energy LFV decays. In Section 5, we discuss in more detail the cross sections, comparing 
with the expected signals at the LHC, where appropriate. Finally, in Section 6, we summarize 
our conclusions and outlook. 



2 LFV in Slepton Pair Production 

2.1 CMSSM with LFV 

In the unrotated charged- lept on flavour basis £i = (e^, /11, tl, e* R , p,* R , f R ), the charged slepton 
mass matrix is: 

/ Mh M? R \ 

where 

M 2 LL = mlm e + M 2 L - l -{2m 2 w -m 2 z ) cos 2(3 I, 

M RR = m\m t + M R - (m| - m 2 w ) cos 2(3 /, 
m lr = ( Ae ~ A* tan (3) m e , 

m rl = (M 2 LR )l (2) 

Here we parametrize trilinear soft supersymmetry-breaking terms as A|-(A e )jj, where the A e are 
the respective Yukawa couplings. For universal soft terms at some high input scale, one has 

M\ = M 2 R = ml I, A%=A^ l3 , (3) 

whereas flavour-mixing entries may be parametrized by: 

4x = {M 2 xx f/{M 2 xx f (X = L,R). (4) 

The correspondence between the mixing parameters 5xx and the flavour mixing parameters 
used in the phenomenological study of slepton production cross sections is easily derived. For 
example, defining the splitting of the third-generation soft supersymmetry-breaking mass terms 
as might be generated by the renormalizat ion-group equations (RGEs) by 

K XX = 1- (M 2 xx r/(M 2 XX ) 22 , (5) 

the slepton mixing angle #23 and the splitting of the second- and third-generation mass eigen- 
states Am 2 3 can be written as: 

tan(2^ 23 ) = 25 XX /K XX , (6) 

Am 23 = (7) 
m sin(2fl 23 ) 

3 





Figure 1: Feynman diagrams for e + e — > — > ^t^ a x\x\- The arrows on scalar lines indicate the 
flow of lepton number. Similar diagrams - appropriately modified - exist for charginos. 



where m = | (777,2 + ^3)- 

The evaluation of the LFV observables is done by performing the diagonalization of the slepton 
mass matrices (see [8], for instance), inserting the full rotation matrices in the lepton-slepton- 
gaugino vertices and summing over all the mass eigenstates of the exchanged particles. 



2.2 LFV Cross Sections 

As already discussed in the introduction, charged-lepton flavour violation at a LC may oc- 
cur either directly in slepton pair production or indirectly via slepton production in cascade 
decays [TO] . Processes leading to lepton production in the decays of a pair of sleptons include 

e+ e - _> Irlf^r+^xlxl 

e + e- ->• vo)] ->• ^n^xtxl, (8) 

for which representative Feynman diagrams are shown in Fig. 1. Slepton production may also 
result from the the cascade decays of the heavier gauginos, e.g., via the processes 

e + e~ -> ^^r^^x?. (9) 
These amplitudes differ in the interference terms, and read as follows: 

for pair production and cascade slepton production, respectively, where Mp and Md are the 
production and decay amplitudes for sleptons in the absence of LFV, and T ia parametrizes the 
lepton- flavour- mixing matrix element. 
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Complete expressions for the respective cross sections are given in [TT], and used in our work. 
However, intuition may be developed by first looking at the simplified case including only 2-3 
slepton mixing. In this case, the only diagrams contributing to the cross sections are those 
mediated by 7 and Z° exchange. Even simpler formulae can be obtained in certain limits, 
such as the cases of narrow widths and small mass differences between the sleptons of different 
generations. For example, if Afhij <Cm = \{jn 2 + m 3 ) and mT^ ~ {fhjTi + rhjTj)/2 <C rh 2 , in 
the presence of 2-3 mixing the cross sections for pair and cascade slepton production are given 
by [22]: 

<x pair = x 23 (3 - 4 X23 ) sin 2 2^ a(f f -> t t)Br(£ + -> t X)Br{£~ -+fY), (12) 
a c a f = X23 sin 2 2fif 23 a(f f -> £ + X t)Br{t l~ Y\ (13) 

where X,Y represent gauginos in the final state, and cr(f f — > £ + £~) and _Br(£ ± — > £^X) 
are the cross sections and branching ratios in the absence of flavour violation. In the case 
of flavour conserving cascade decays, a(f f —> £ + X £~) includes the intermediate production 
of two gauginos and the subsequent decay of the heavier into a lepton-slepton pair. LFV 
enters via X23 = ^23/2(1 + ^23) where x 2 3 = Am 2 3/r 2 3, and sin 2 2023 parametrizes the slepton 
mixing angle. In the limit X23 S> 1, interference can be neglected and the results are even 
further simplified. However, in the case that the interference term dominates, it suppresses 
LFV significantly. 

The exact experimental signatures of the above effects depend on the decay chains, as discussed 
in previous analyses [IT] . Channels involving the lightest charginos are potentially interesting, 
particularly those with hadronic chargino decays, since the lightest neutralino is 'seen' as missing 
energy. In general, the dominant chargino decay is xt ~^ whereas channels such as xl ~^ 
QuQdX are relevant when < or, in other words, when the xt is the NLSP. Channels 
leading to hadronic xt decays [TTJ are mediated only by left sleptons, whilst e + e~ — > r ± /i T +2x° 
may be mediated by either right or left sleptons. However, points with m^+ < are disfavored 
in the framework of the CMSSM. For instance, the previously proposed point with tan — 3, 
m = 100 GeV, M 1/2 = 200 GeV [J2] was excluded by the LEP m h bound [18] and is now 
incompatible also with LHC limits [19]. This could perhaps be achieved in the focus-point 
region, where M\n < 1 TeV and mo > several TeV, but slepton masses in this region are very 
heavy and hence LFV processes are more suppressed. 

The channel e + e~ — > r ± /i =F + 2x° was studied in [T3], where it was found that a signal of ~ 1 fb 
could be distinguished from the background. Only slepton mixing in the left sector arising from 
see-saw neutrino mass models was included in |14j . and the limit from r — > /ry was at the time 
less restrictive. As shown in [17], reasonable rates of LFV processes at the LHC in the slepton 
sector are possible while preserving the bound on r — > fi'y if one assumes mixing also among 
the right sleptons, as may occur in non-minimal GUT theories. 
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Figure 2: Contours of the cross sections a(e + e~ — > l + l~) = 1 fb in the (M]y2> m o) planes for 
tan/3 = 35 (left) and 45 (right) at three different energies, namely y/s = 500 GeV (red thin lines), 
y/s = 1000 GeV (black thick lines) and y/s = 2000 GeV (blue thin lines). The upper and lower 
near-horizontal lines correspond to the following slepton pair production processes, respectively: t{~ 
(solid) and v T v T (dash). From left to right, the near-vertical lines correspond to a{e + e~ — > X1X2) = 1 f° 
(dot-dash) and a(e + e~ — > X1X2) = 1 ft ( double- dot- dash) . In regions above the line where m~+ = mp 
(purple dots), chargino-mediated processes are allowed. Representative contours of BR(b s^) are 
indicated, and the regions with LSP relic density in the range preferred by WMAP are shaded green. 
The benchmark points selected for further study are indicated by crosses. 

3 Supersymmetric Parameter Space 

From the above discussion, it is clear that the dependence on the supersymmetric parameter 
space is particularly crucial, since the chains of cascade decays that will dominate, as well as 
the relative strengths of the signals and backgrounds, depend on the mass hierarchies of the 
superpartners. 

We start with some preliminary considerations of the relevant channels within the CMSSM pa- 
rameter space. In previous work [T7J [23] a region that is promising from the cosmological point 
of view was discussed, in which sizable f — \ coannihilations or direct-channel H/A resonances 
lead to values of Q x h 2 consistent with WMAP. As a representative of the coannihilation case, 
we chose tan = 35, a value that leads to a relatively heavy sparticle spectrum compatible 
with WMAP. At tan f3 = 45 we found points in the resonance funnel, which have larger scalar 
masses since mo can take higher values than in the coannihilation area. We did not consider 
the focus-point region of the CMSSM because the super-heavy slepton masses there lead to 
reduced LFV effects. We note that these regions are compatible with the LHC constraints on 
supersymmetry, which indicate that Mx/% > 450 GeV and favour tan > 15 [20] . 

The physical values of the masses are obtained by integrating the renormalization group equa- 
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tions from Mqut down to low energies. We use the two-loop supersymmetric renormalization 
group equations [24] except for the trilinear terms, the gaugino and the sfermion masses, which 
are calculated at the one-loop level. Electroweak Symmetry breaking (EWSB) occurs at the 
scale Msusy = y/m tl rh t2 (where the one-loop contributions to the scalar potential are less 
relevant). At this scale we incorporate the SUSY threshold corrections to my, , m T and mt by 
redefining the corresponding Yukawa couplings as done in ref. [25J. From Msusy to Mz we 
use the Standard Model renormalisation group equations. The running top mass is calculated 
iteratively, removing the Yukawa coupling At, along with its derivative, from the remaining 
running from m t to Mz- 

Within this framework, for every set of input parameters, we perform an iterative integration 
of the RGE's (back and forth, from high to low energies). We define Mqvt as the meeting 
point of the gauge couplings «i and ot 2 and use this scale to find the unified value of the 
couplings, otjj . The value a 3 (M GUT ) is obtained from a 3 (M Mz ) (its deviation from ay being 
very small). The third generation Yukawa couplings are obtained by using the following set 
of values: m T {Mz) = 1-7463 GeV (which takes into account the SM radiative corrections) 
and the top pole mass m t = 172.6 GeV [26]. For nib we take nibiMz) = 2.92 GeV which 
is the value where m(m,b) = 4.25 GeV is mapped, taking a s (MZ) = 0.1172. Note that the 
value of the bottom Yukawa coupling extracted from nibiMz) is very sensitive to the SUSY 
threshold corrections, especially at large tan/3; in our case we implement a complete one loop 
contributions as in ref. |27j. 

In Fig. [2] we present contours of the production cross sections for sparticles that are relevant 
for our study of LFV effects within the CMSSM , using 1 fb as a reference cross-section value. 
We indicate in green the constraint on the parameter space derived from the WMAP-favoured 
CMSSM parameter space (0.094 < Vth 2 < 0.128 at 3 - a) using MicrOMEGAS [28], and in grey 
the regions excluded by the constraints that BR(6 — > 57) > 2.5 • 10~ 4 and nih > 111 GeV ( we 
assign an uncertainty of ~3 GeV to the theoretical calculation of m/,). We present contours for 
the production of r+ff" calculated using CalcHEP [25]: the cross sections for producing other 
right sleptons would have similar contours, though at somewhat smaller values of (Mi/2,mo) 
because they are heavier than the f\. The left slepton production cross sections are smaller, 
in general, and here we draw the o{y T v T } = 1 fb contour (the cross sections for all slepton 
pair-production processes are compared in Fig. [3]). The relevant channels for LFV in cascade 
decays following pair production are also displayed in Fig. [21 along with the line m-+ = m^, 
above which chargino-mediated LFV processes are allowed. 

We would also like to note the following, regarding comparisons with known results: For tan = 
30 our results compare very well with the ones obtained using the SUSPECT code [30J, as 
implemented in Micromegas |28] . At tan/3 = 45 we observe very small differences, due to the 
fact that we include all SUSY threshold corrections to nib and m T (from ref. [25]). Note that 
the resonances on the \ annihilation channels, relevant at large tan/3, are dependent on the 
value of Xb, and therefore on the accuracy of the approximation in computing nib- To quantify 
the difference, we note that our results for tan/3 = 45 are similar to the ones obtained with 
suspect for tan/3 = 46 and m 6 (m b )=4.4 GeV. 
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Focusing initially on the lines corresponding to yfs = 500 GeV, we see that on-shell production 
of left sleptons would occur mostly in the region excluded by 6 — > 57. At this energy, therefore, 
we expect right slepton production to be more important. However, production of charged 
slepton pairs in the cascade decays of neutralinos is also possible, whilst chargino decays are 
disfavored. At y/s = 1000 GeV, we see a significant part of the parameter space with left slepton 
pair production (below the dashed line) and with both left and right slepton production (below 
the solid line). The production of pairs of charged sleptons in chargino cascade decay is not 
possible in most of the WMAP favored area, since m-+ < m^. Finally, at = 2000 GeV all 
the decay channels are present in the region of parameter space favoured by WMAP. 

Based on this discussion, two benchmark points in the allowed regions for each of the choices 
tan/3 = 35,45 have been selected for subsequent study, namely: 



(a35) 


tan /3 


= 35 


m 


= 255 GeV, 


M 1/2 


= 610 GeV, 


(635) 


tan j3 


= 35 


m 


= 345 GeV, 


Mi /2 


= 900 GeV, 


(a45) 


tan j3 


= 45 


m 


= 440 GeV, 


M 1/2 


= 550 GeV, 


(645) 


tan/3 


= 45 


m 


= 730 GeV, 


M 1/2 


= 900 GeV 



with Aq = in each case, resulting in the particle spectra shown in Table 1. We note that 
all these points lie beyond the regions of CMSSM parameter space excluded by 2010 LHC 
data [TjJJ [2D]. In general, these points lead to more promising detection prospects at a LC at 
high energies: at a/s = 500 GeV, only production through right sleptons is allowed. 



Point 


m 


Ml/2 


tan /3 


M 4 






M xl 


M fl 


M f2 


M ln 


m !l 


My r 


Mv e 


a35 


255 


610 


35 


480 


741 


255 


480 


261 


477 


342 


478 


454 


472 


b35 


345 


900 


35 


720 


1034 


383 


720 


385 


670 


478 


683 


655 


679 


a45 


440 


550 


45 


432 


671 


229 


433 


338 


543 


485 


571 


519 


566 


b45 


730 


900 


45 


725 


1015 


385 


725 


590 


872 


801 


936 


858 


933 



Table 1: Relevant sparticle masses (in GeV) for the selected benchmark points. The selected points 
all have Aq = 0. 



We show in Fig. [3] the dependence of the flavour-conserving (FC) cross sections on the \fs of the 
collision for the benchmark points ffl4"]) with tan /3 = 35. The cross sections for the benchmarks 
with tan (3 = 45 are very similar. We can already get an idea of the LFV cross sections from 
the magnitudes of these FC cross sections, by looking at the approximate formulas ( fl2|) and 
( fl~3i) . which assume that only 2-3 mixing is present. These provide an initial picture of the 
order of magnitude of LFV (up to 40% error), simply by convoluting the slepton flavour- 
violating decays with the FC cross sections. The latter can be extracted automatically, e.g., 
from CalcHEP [29], allowing one to explore quickly the relevant parameter space. If the mixing 
was indeed limited to the 2-3 sector, the accuracy would be larger, within 10% of the exact 
calculation discussed at a later stage. 

As shown below, taking into account the full flavor mixing in the vertex [HI [11] leads to additional 
channels and an enhancement in the cross sections a(e + e~ — > iilj), also for the case of mixing 
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Figure 3: Energy dependence of the cross sections for sparticle production in the channels (Ejj and 
(GJ). The black thick (blue thin) lines correspond to point a35 (b35) in ji^p . The upper left (right) 
panels display left (right) charged slepton production, and the lower left (right) panels display sneutrino 
(gaugino) production. 

in the LL sector, despite the fact that the mixing with the first generation is very constrained 
by the recent MEG bound BR(/i -> e^)< 2.4 • 10" 12 [31]. 

The cross sections e + e~ — > £~£f, i — L, 1, R, 2, which are relevant for calculating e + e~ — > 
i^ij' — > T ± /i =F XiXi) are displayed on the upper panels of Fig. [3j We see that the production 
rates for left and right sleptons are comparable. 

In the case of e + e~ — > v(i) c - — > ^^XiXi^ Fig. |3] indicates that the production rate for sneu- 
trinos can be larger than that of charged sleptons. However, i) even if e + e~ — > Vii) c - is large, 
the decay z>j — > is allowed only in a restricted area of the parameter space, as we can see 
in Fig. El and ii) the advantage of producing the Xi is n °t so important in this case, because 
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the decay Xi ~ ► X + QuQd is not allowed in the range of tan (3 we favour (we recall that in the 
CMSSM the NLSP is usually a stau). 

Finally, we comment on the cascade slepton production processes of eq. (Q. These cases are 
not so promising within the CMSSM, but the analysis is simple, using the flavour- conserving 
cross sections for charginos and neutralinos and multiplying by the LFV branching ratio. The 
LFV signal arises from the decay of the X2 or X2 > which are produced at lower rates than the 
sleptons. The LFV decay of the X2 has been already studied in the context of the LHC in [17] . 
Here, in agreement with [TT] . we find that mixing in the RR sector may lead to larger signals. 
However, the resulting cross sections are smaller than those mediated by sleptons. Also, the 
channel with xt is not so relevant, not only because the production of charginos is lower than 
that of neutralinos (as can be seen in Figs. |3]), but also because it is only sensitive to LL mixing 
in points where m~+ < . 

Other production channels arising from xfx^ and XiXji where i, j denote heavier neutralinos, 
are also open at higher energies and can contribute to the LFV signals. However, at the points 
we consider their production is lower than the one we study in detail and do not give any 
additional insight. 

4 Opportunities for LFV observation at the LC 

LFV effects are severely constrained by decays of the form BR(£j — > £jj), which limit signifi- 
cantly the allowed parameter space, with respect to both the mixing terms 5ll,rr and the values 
of the soft supersymmetry-breaking masses. These limit the range of parameters in which we 
can get observable cross sections without violating the LFV bounds. In general, we may expect 
that LFV decays are progressively suppressed for a heavier sparticle spectrum at fixed values of 
8ll,rr, and hence that one may get observable effects only for larger values of 8ll,rr- However, 
since the LC gives direct access to energies that are relatively high, the suppression induced by 
a heavier spectrum may be smaller than those in rare LFV decays and /1 — e conversion. It is 
natural, therefore, to ask the following questions, (i) For a heavy sparticle spectrum, how small 
are the 8ll,rr that are accessible before the cross sections become too small to be observable? 
(ii) For a light sparticle spectrum, what are the minimum values of Sll^rr that can be probed? 

These issues are addressed in Fig. HJ where we display values of the flavour- violating parameters 
8rr and 5ll-, showing lines of constant a(e + e~ — > r ± /i T + 2x°) for y/s = 2000 GeV0. Using 
the bounds on r /i7 and r — > /ry from [32], we display the allowed ranges in ((<5ll)23, (8ll)i3) 
planes (upper panels), and the corresponding mixing parameters in the RR sector (lower pan- 
els). In this exploration, the values of (5ll,rr)\2 were chosen to be proportional to those of 
(<*> ll,rr) 13, so as to reduce the number of parameters. We have checked that the dependence 
of our results on these parameters are small, provided they are taken in a range such that the 

1 We choose ^fs = 2000 GeV in Fig. [4] because, had we chosen y/s = 1000 GeV, only the RR channels would 
lead to on-shell sleptons. Even for y/s — 2000 GeV, we can see in Fig. [5] that point 645 predicts very low cross 
sections for the LL case, due to the fact that the left sleptons are only marginally accessible, kinematically. 
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( 8 Rr) 2 2 



Figure 4: Constraints on the magnitudes of the mixing parameters and possible LFV effects for points 
with tan/3 = 35 (left panels) and tan/3 = 45 (right panels). The shaded areas are those allowed by 
current limits on BR(t — > cy) (dot-dash line) and BR(t — > fij) (dash line) using a35, a45 as reference 
points (thick lines bounding the solid shaded areas) and , 635, 645 (thin blue lines bounding the ruled 
shaded areas). The solid lines are contours of o~(e + e~ — > r /x T + 2x°) in fb for ^/s = 2000 GeV. 



bound on BR(/i — > ej) is preserved. In addition, we note that the ratio (5ll,_rr)i3(5ll,_rr)i2) 
has a very small influence in r — > (r — > 

As we see in Fig. HJ the constraints from the LFV branching ratios result in horizontal and 
vertical lines in the 5ll and 5rr planes, corresponding to a box of allowed parameters that 
scales with the masses in the spectrum. On the other hand, the lines of constant LFV cross 
sections depend strongly on the kinematical properties of the chosen point as well as on the 
size of the mixing. In particular, we note the following: 

• The shapes of the lines are quite different for each benchmark point, as was to be expected 
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Figure 5: Values of a(e + e~ — > t^t- — > t^il^ + 2^°) ws /or t/ie benchmark points with tan /3 = 35 
(7e/tj and tan/3 = 45 (right): solid (dashed) lines correspond to a35 and a45 f 635 and 645J. Along 
the thick lines we assume the RR mixing of U6\) , and along the thin lines the LL mixing of i!5\). 

from the variations of the cross sections between the different benchmark points in the WMAP 
favoured area. For instance, for LL mixing, the differences between the predictions for the two 
points at tan j3 = 45 (upper right panel in Fig. Hj) can be understood from the behaviours of 
the cross sections shown in Fig. [51 

• Except for LL mixing at tan/3 = 45, we obtain LFV cross sections in the fb range at points 
of the parameter space where LFV decay branching ratios are below the present experimental 
bounds. This reflects the observation that passing to a heavier spectrum increases significantly 
the allowed size of the mixing parameters, due to the rapid decrease of the LFV branching 
ratios for heavy sparticle spectra. 

• It is interesting to remark that at the LC mixing in the LL sector is observable without 
exceeding the bounds on £j — > £^7 J3TJ [32], due to the presence of additional channels as 
compared to the LHC In the case of the LHC [17), mixing in the RR sector (which is less 
constrained but would require a departure from the CMSSM) had turned out to be the most 
promising avenue for discovering flavour violation. This is a clear distinction between the 
two colliders, and implies an advantage for the LC, in that its searches can be more directly 
connected with neutrino mass and mixing parameters. 

The above results motivate a particular choice of parameters for a more detailed discussion of 
LFV cross sections. We assume the same values of the mixing parameters for both choices of 
tan/3, namely: 

(Sll)i3 = 0.02, (5 LL ) 23 = 0.02, (15) 
OMia = 0.04, (<W>23 = 0.15. (16) 
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For completeness, we also introduce a small mixing between the first and second generation: 
($ll/rr) 12 = 0.2 ■ (5ll/rr)vz- The LFV decays into r — \i pairs are not heavily dependent on 
this parameter while, with this choice, the bound on BR(/i — » ej) does not over-constrain the 
parameter space. 

LFV in the \i — e channel is generally suppressed in our scheme, unless Sl R mixing is introduced. 
In our analysis, however, we do not include such a mixing, in order to be able to discuss 
separately the LFV arising from each chiral sector, thereby avoiding the additional model 
dependence associated with choosing 5lr 7^ 0. In the case of large RR mixing, non-trivial 
interference of the RR sector with the LL sector would arise through LR mixing. The \i — > 
bound, which could almost be ignored in selecting RR and LL mixing, would now play a very 
important role, forcing us to either a very small part of the parameter space or to very small 
values of (Sl R ). To reach the bound on radiative decays with LR contributions alone we would 
need: (5 LR ) 12 ~ 10~ 3 , (6 LR ) 13 ~ 10" 2 [33]. 

We also followed the conservative approach of keeping A = 0, avoiding the introduction of 
other new LFV parameters. The main consequence for our analysis of allowing Aq 7^ would 
be to allow the chiral composition of the staus to vary, since the first and second generations 
have small Yukawa couplings. However, additional LFV could be introduced by assuming a non- 
trivial mixing matrix for the A terms. In this case, chiral mixing would occur for all generations, 
restricting among other effects the size of the RR mixing. The decay fi — > cy is the most 
sensitive to LFV A terms, and establishing a correlation between \x and r decays without having 
a specific model would be non-trivial. Nevertheless, since the bound on BR(/i — > ej) imposes a 
very severe restriction on (M LR )i 2 , even assumptions such as (M LR ) 23 ~ (M LR )i 3 ~ 100(M M ) 12 
would not alter our main conclusions. 



5 Study of LFV at the LC 

The next step is to understand the dependence of the results on the available y/s. Clearly, at 
the LC, higher values of y/s provide accessible cross sections for heavier sparticle spectra, a 
potential advantage of the LC over the LHC, where it is relatively difficult to reach large values 
of the scalar masses. The variation of a(e + e~ — > r ± /i =F + 2%°) with y/s is displayed in Fig. O 
where it is shown how, as ^fs = 2P cm increases, it becomes possible to produce on-shell sleptons 
with larger masses. We should underline once more that the effects due to the kinematics are as 
relevant as those from the LFV mixing parameters. The shapes of the lines can be understood 
as follows. 

• The cross sections for slepton pair production show interesting dependence on y/s, since addi- 
tional channels arise from slepton mixing. We see in Fig. [5] steps in the left slepton production 
cross sections, which appear as the sleptons of the first two generations become kinematically 
accessible. The steps corresponding to lighter generations are less pronounced in the case of 
RR mixing than in the LL case. 

• The overall cross sections for slepton production in cascade decays are independent of the 
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Figure 6: Same as in Fig. but for cascade slepton production in the channel a(e + e — > X1X2 ~* 
t ± /^ + 2 X °). 

LFV parameters, and LFV arises in the decays of the gauginos. Figs. |6] and [7J display the cross 
sections for the pair production of X2X1 or X2X1 multiplied by the LFV branching ratios for 
the decays of x° or X21 computed as in [17]. This explains the simple shapes in this case. 

• Distinguishing between pair and cascade slepton production is quite crucial, particularly since, 
as already mentioned, it is a fundamental difference between the LC and the LHC Comparing 
the pair and cascade slepton production cross sections in Figs. |5] and |6] (where we chose similar 
scales to facilitate comparisons) we see the following: in slepton pair production, both LL 
and RR mixing effects are of the same order of magnitude; on the contrary, in the case of 
cascade slepton production, the channel with LL mixing has a cross section at least one order 
of magnitude less than that with RR mixing, offering a distinction between the two mixing 
scenarios. 

• Channels with charginos in the final states could be promising when < m^, i.e., when 
the xt is the NLSP. However, as already remarked, this is not the case of the CMSSM. In 
general, the dominant chargino decay is xt ~^ thus we do not obtain additional information 
beyond that from simple slepton production and decay. This is true also for channels with the 
heavier neutralinos. 

• Cascade slepton production mediated by charginos is possible only with mixing in the LL 
sector. Fig. [7J (left panel) indicates that this channel could be relevant at ^fs larger than in 
the neutralino mediated channels, and the cross sections would be of comparable magnitude 
in the LL mixing case. On-shell production requires 777.5 > ^v+> which is not the case for the 
sample points with tan p = 35, as they lie in the coannihilation region of the WMAP allowed 
area. However, for tan = 45 the WMAP region includes points with resonant annihilation of 
neutralinos where the condition mr, > m~+ is satisfied. 
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Figure 7: Left panel: Same as in Fig. but for cascade slepton production in the channel a(e + e~ 
xtx} ~* t^I^ + 2x°). Right panel: Same as in Fig. but for pair production with sneutrino 
mediation o~{e + e~ — > viu^ — > r =t ^ =F + XiXi)- F° r tan/3 = 35 these channels are not allowed because 
mo < . 

• The expectations for sneutrino pair production and subsequent chargino production are shown 
in the right panel of Fig. [7J Relaxing the assumptions of the CMSSM can change the kinematics 
of e + e~ — >■ XiX^i as we U as the mass ratio of xt versus v that determines the decay xf ~^ 
Xi = T ± yU =F . The branching ratios of the xf decays to u T fi are 80% and 65% at points a45 and 
645 respectively, with a subsequent ~ 100 % decay of the stau into a r and a Xv 

Cross sections for e + e~ — > r ± yU =F + 2x° above 1 fb can be reached in the areas of parameters 
shown in Fig. [8] for energies below those indicated. The large areas between the thick solid 
red, black and blue lines, on the one hand, and the shaded regions that are excluded by the 
indicated present experimental constraints, on the other hand, demonstrate that there are 
ample opportunities for LFV discovery and measurement at the LC. These opportunities are 
exemplified by the benchmark points marked by crosses in Fig. [HJ 

In Fig. [9] we present the maximum values of the cross sections in the areas allowed by all the 
current constraints. The shaded areas show the possible ranges of <r(e + e - —> + 2%°), 

and we note that points along the WMAP strips (solid lines) generally have high values of the 
cross sections. The dashed lines show the possible values of cx(e + e~ — > r ± e =F + 2%°) along the 
WMAP strips, and we see that these cross sections may be of the same order of magnitude as 
those for /i — r pairs 0. 

In Fig. [10] we display the expectations for the production of LFV r — /i pairs in the cascade 
decays of eq. (EJ), which are suppressed compared to the pair production process of eq. (jHJ). 

2 We do not display values for a(e + e~ — > /j + 2%°), because this would require a more precise analysis of 
mixings in the LR sector. 
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Figure 8: The solid red (black) (blue) lines are contours where a(e + e~ — > i ^ + 2%°) = 1 fb at 
y/s = 500, 1000, 2000 GeV, assuming the left mixing parameters U5\) in the upper panels and the 
right mixing U6\) in the lower panels. The left panels are for for tan (3 = 35 and the right panels 
for tan/3 = 45, assuming A$ = in both cases. Each panel also shows the areas excluded by current 
bounds on BR(/i — > ej) (thin solid line), BR(t — > e^f) (thin dash line), BR(t — > ^) (thin dot-dash 
line). The areas excluded by BR(b — > sj) and the LEP Higgs search are also displayed, and the green 
area denotes the WMAP favored region. We see that there are ample opportunities for LFV discovery 
and measurement at the LC, and the benchmark points chosen for further studies are indicated by 
crosses. 

The production of LFV pairs from x° decays is analogous to the corresponding process at the 
LHC [T7]. We observe that phenomenologically interesting values on the LL sector induce very 
low cross sections, whereas mixing in the RR sector may be more observable, as it enhances 
the expected rates. 

The SM background can be as large as several tens of fb. However, it has been shown in ref. [H] 
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Figure 9: The shaded areas show the possible ranges of a(e + e~ — > t ^ + 2%°) for tuq < 1000 GeV, 
tan/3 = 35 (left panels) and tan/3 = 45 (right panels), assuming Aq = in both cases, with ^fs 
fixed to 500 GeV (grey), 1000 GeV (green) and 2000 GeV (orange). In the upper panels, we assume 
the left mixing parameters ( fi5j) . whereas in the lower panels we assume the right mixing (16]) . The 
solid lines present the possible values for models along the center of the WMAP strips, with the lines 
corresponding to ^/s=500, 1000, 2000 GeV being progressively thicker. The corresponding predictions 
for a(e + e~ — > r^e^ + 2%°) are shown by the dashed lines. 

that after imposing several cuts, it is possible to distinguish LFV signals at the level of 1 fb. 

The SUSY background arises from flavor conserving channels like e + e~ — > ViOi^ e + e~ — > 
and subsequent cascade decays that result in r ± /i T + . This background is very dependent 
on the SUSY spectrum and is much smaller than the SM one. In our case, the branching ratios 
for e ± r T + |£and /i ± r T + arising from the production of e£e^ and ft^^L are almost zero at 
the points with tan/3 = 35, and of about 10% at the points a45 and 645. Without considering 
any cut, this induces a background of less than 1 fb in r ± /i T + ^at energies below 2 TeV. At 2 
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Figure 10: Predictions from cascade decays (GJ) along the WMAP strips for yfs= 1000 GeV (thick 
black lines) and 2000 GeV (thin blue lines). These channels are not accessible at ^/s=500 GeV. The 
solid (dashed) lines correspond to mixing in the RR (LL) sectors. 

TeV, the background for e ± r =F + ft is bigger due to the higher production of and v in e ± e T 
collisions (while at 645 it remains of the order of 1 fb, it reaches a value of 12 fb at the point 
a45). Other channels like e + e~ — >■ xtxj an d e + e~ — >■ xt e ~^e produce backgrounds at the fb 
level even at 2 TeV. 



6 Conclusions 



Motivated by the challenging possibility of observable charged lepton flavour violation, which 
would provide a new window on flavour physics that illuminates novel aspects of supersymmetry, 
we have explored possible signatures at the Linear Collider (LC), combining phenomenological 
constraints with the cosmological considerations implied by WMAP. We have found that both 
direct slepton pair production and sleptons produced in cascade decays may provide interesting 
signals in the cosmologically-favoured region of the supersymmetric parametric space. More- 
over, the LC could provide additional insights beyond those obtainable from the LHC, by virtue 
of its greater kinematic range for slepton production and its sensitivity to LL mixing as well 
as RR mixing. 

Within this framework, we found the following. 

• The LC enhances significantly the prospects of detecting LFV for heavy sparticle spectra, 
where flavour-violating rare decays and conversions are significantly suppressed. This allows 
probing an entirely different range of the flavour-violating parameters. 

• Unlike the LHC of ref. [17], at the LC we can have significant LFV within the CMSSM via 
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LL mixing. 

• The decay e + e~ — > r ± /i T + 2%° provides an optimal search channel in this respect. Channels 
where charginos play a significant role may require a departure towards theories where gaugino 
or Higgs unification is broken, in order to respect the conditions imposed by cosmology. 

• Comparing direct slepton pair production with indirect slepton production through the cas- 
cade decays of heavier sparticles, we find that in the latter case, the cross sections induced by 
LL mixing are approximately one order of magnitude lower than in the RR mixing case. This 
is consistent with the results previously obtained for the LHC in Ref. jl7j, and is due to the 
presence of additional channels in cascade pair production. These lead to an interesting scaling 
of the cross sections with the available energy, which is sensitive to slepton mixing parameters 
of the three generations. 

It would be very interesting to further investigate the following issues. 

• The cross sections expected in non-minimal extensions of the theory might not only enhance 
channels that in the current scheme are more suppressed, but could also enable a comparison 
of the allowed range of mixing parameters in different models. 

• The expectations for models of massive neutrinos, in which quantum corrections provide a 
significant source of LFV in the LL channel, provide a potential link between LC observables 
and neutrino mass and mixing parameters. 

• Detailed simulations of signals and backgrounds are desirable. 

Overall, it seems that the LC provides an optimal environment for the study of LFV, whereas 
the LHC is limited to specific channels that have significant backgrounds. The fact that the LC 
opens up additional possibilities may prove significant for making the link between observable 
cross sections and flavour model building. 
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